Enlightening complexity: making energy with chaos 
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We study the energy harvesting of pliotons undergoing chaotic dynamics with different com- 
plexity degrees. Our theory employs a multiscale analysis, which combines Hamiltonian billiards, 
time-dependent coupled mode theory and ab-initio simulations. In analogy to classical thermody- 
namics, where the presence of microscopic chaos leads to a single direction for time and entropy, an 
increased complexity in the motion of photons yields to a monotonic accumulation of energy, which 
dramatically grows thanks to a constructive mechanism of energy buildup. This result could lead 
to the realization of novel complexity-driven, energy harvesting architectures. 
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Introduction. — One of the most important theo- 
retical dispute of the 19th century concerned the foun- 
dations of statistical mechanics and was pointed out by 
Loschmidt in a famous paradox [1|. He objected, in par- 
ticular, that it was not possible to deduce an irreversible 
process from the time-reversible dynamics provided by 
Boltzmann theory, which was based on the motion of 
classical particles undergoing elastic collisions . As ex- 
perimentally ascertained in modern times ^ , the physi- 
cal machinery behind this apparent contraddiction is the 
presence of microscopic chaos, which exhibits an irre- 
versible behavior due to the unpredictability of the re- 
sulting dynamics. Sensitivity to the initial conditions re- 
lating the thermodynamic system and its environment to 
the microstate level, in fact, gets exponentially amplified 
by chaos and culminates into completely irreversible tra- 
jectories despite the time reversibility of the microscopic 
Hamiltonian. The consequence of this result is notewor- 
thy and remarkable, as it establishes a single direction for 
both time and entropy [4] . This is not only at the basis of 
our everyday perception of the world, but has also signifi- 
cant implications in several fields of research ranging from 
quantum mechanics to cosmology [1, @] . In the language 
of complex systems, the existence of chaos, or unpre- 
dictability, is the signature of a nonvanishing complexity 
in the dynamics. Complexity stems from a nonzero en- 
tropy and may be quantified by the Kolmogorov-Sinai 
entropy of the corresponding chaotic motion ^] . In light 
of this interpretation, we may conclude that the influence 
of a varying complexity even on a simple dynamics, such 
as the one considered by Boltzmann, may yield dramatic 
consequences on the resulting macroscopic observables. 
On the other hand, the classical formulation of statisti- 
cal mechanics was originally conceived to study ensem- 
bles of atoms and molecules interacting through electron 
wavefunctions, and is therefore characterized by the oc- 
currence of short-range repulsion among the system con- 
stituents [8]. With these premises, a natural question 



therefore arises and concerns the effects of complexity 
on particles ensembles showing a bosonic behavior, such 
as photons, whose nature does not entail any repulsion. 
Contrary to electrons, photons are energy-carrying parti- 
cles (of multiples of the light quantum hoj) and their fun- 
damental macroscopic observable is the energy. In this 
respect, our initial question becomes: what is the energy 
harvesting potential of a complex system of photons ? It 
is not difficult to realize the significant importance of this 
problem. The theme of energy harvesting, in fact, is one 
of the leading challenges of contemporary science 
and possesses crucial implications which go beyond ap- 
plied physics and largely permeate fundamental research 
on living organisms 3, 3|- 

In this Letter, we address the above mentioned open 
question. We begin by considering the simplest inte- 
grable, time-reversible Hamiltonian modeling the clas- 
sical dynamics of photons and then add integrability- 
breaking terms providing complexity. We then analyze 
the energy harvesting capacity of the system when pho- 
tons are injected from an external environment. Al- 
though there exist a large literature on the relaxation 
of energy inside resonators [l8| . the problem of energy 
harvesting has never been challenged. The problem is 
tackled by employing a multiscale analysis. At variance 
with "simple" systems, in fact, complex dynamics shows 
the contemporary presence of different scales, each char- 
acterized by a specific set of fundamental laws [15| . A sin- 
gle scale highlights just a facet of the phenomenon, and 
all scales are required to gain a complete understand- 
ing over the problem. With reference to the propaga- 
tion of photons, we identify three main spatio-temporal 
scales: i) the wave-carrier dynamics, which yields the 
most detailed description and it is based on the full set 
of Maxwell's equations; ii) the wave-envelope time-scale, 
which provides a dynamical description at a larger time 
with respect to photon's inverse frequency iii) the 
short- wavelength limit, which constitutes the most ab- 
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FIG. 1. (Color Online), (a) Complexity C versus symmetry 
parameter a; (b)-(c) open billiard geometry for a = (b) and 
a = 0.5 (c). 



stract scale and entails a classical description as embod- 
ied by Hamilton's equation [Tij , applied when both time 
and space are large with respect photon wavelength A 
and characteristic time respectively. We begin with 
the short- wavelength limit, where the problem is math- 
ematically defined, and then analyze the dynamics on 
finer scales up to the most detailed level. The main re- 
sult of this work is the demonstration of a cooperative 
dynamics among photons, originating from the complex- 
ity of their classical motion and sustaining a construc- 
tive mechanism of energy buildup. The efficiency of this 
process increases with the complexity and, quite remark- 
ably, with the number of photon frequencies injected from 
the environment, thereby opening the way to the realiza- 
tion of broadband, complexity-driven energy harvesting 
schemes. 

Short-wavelength scale: stochastic transition to com- 
plexity. — Stemming from the Maxwell's Hamiltonian 
modeling photons dynamics at scale i): 



(1) 



being E, H electric and magnetic fields, e(r) and /i(r) di- 
electric and magnetic permittivity, respectively, we apply 
the WKB transformation E, H cx e^s{r)-tut^ ^^ms recast- 
ing the original system into the Hamilton- Jacobi (HJ) 
classical form dk/dt = —dD/dr, dr/dt = (9Z?/9k, with 
position X, momentum k = WS and D — Z3(a;,k, r) be- 
ing the dispersion relation of Maxwell's equations. The 
HJ equations model photons dynamics when both spatial 
L and time r scales are large enough, i.e., for |k|L 1 
and UJT ^ 1. Following the ideas of Boffetta et al. p], 
we quantitatively define the phase space complexity C by 
the Kolmogorov-Sinai (KS) entropy: 



A(x, k)(iz^. 



(2) 



with where the standard form of KS entropy is used 
(i.e., Eq. [2]integrated over the Liouville measure), we de- 
fine a normalized version of this quantity by employing 
dv ~ 0(A)/A(ix(ik, with 8 being the Heaviside func- 
tion. The expression obtained is limited in the range 
< C < 1. The complexity parameter C measures the 
relative volume of the phase space (r, k) that encom- 
passes a classical motion with a nonzero entropy. For 
Hamiltonian phase spaces characterized by regular KAM 
regions and chaotic seas, Eq. ([2]) quantifies the growth 
in size of chaotic components and the increase of the un- 
predictability of the overall dynamics. 
In the case of dielectric uniform media, D = £)i(w, k) — 
Lu — ch/ui with Hi being the index of refraction. In this 
condition the HJ equations reduce to that of a free par- 
ticle and the system is integrable. Unpredictability can 
then be supplemented by the use of a symmetry-breaking 
potential. In particular, we use a piecewise uniform di- 
electric geometry by embedding the system into a second 
uniform material Do{uj, h) = uj — ck/uo characterized by 
a different refractive index Uo- The region encompassed 
by Hi defines the harvesting geometry of the problem (i.e., 
where the energy is stored), while the space delimited by 
Uo characterizes the environment interacting with rii. At 
this spatio-temporal scale, the dynamics of photons in- 
jected into the inner region belongs to the class of open 
billiards [17[: the particle represented by (r,k) under- 



goes a free motion and gets reflected at the boundary S 
between and Uo due to the conservation of k; energy 
leaks are then originated from the transmission at iS due 
to Fresnel's laws 18|. In the billiard dynamics, the com- 
plexity of photons is ruled by the symmetry properties 
of S. For our problem, we employed a two dimensional 
polar (p, ip) curve S, with fixed area V and whose shape 
originates from the deformation of a circle: 



p — \/Vl'K — + a cos tA, 



(3) 



being A the largest Lyapunov exponent integrated over 
a measure dv of the phase space. However, at variance 



being a a symmetry-breaking parameter (Fig. [TJd-c). 
Figure [T^ shows the behavior of the complexity C when 
the symmetries of S are changed through a. As a is 
increased from the value a = 0, we observe a marked 
stochastic transition at a w 0.2, with the overall com- 
plexity sharply rising due to the breaking of resonant 
tori. For C < 0.5, the chaotic dynamics is mediated by 
the survival of KAM islands of integrability; at a = 0.5, 
conversely, all tori break and the system shows the max- 
imum degree of unpredictability, with every input condi- 
tion leading to chaos. The numerical evaluation of the 
complexity is carried out by a parallel algorithm, which 
partitions the billiard phase space tp,0 = Zn, k/|k| (n is 
the normal to S) into boxes dipdd = lO^'^rad^, and then 
numerically evaluates the maximum Lyapunov exponent 
A(-0,0). The integration of ([2]) is then performed with 
the trapezoidal method. 

Wave-envelope limit: photons cooperation effects. — 
At the wave-envelope level, we model the system dynam- 
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ics by diagonalizing the Hamiltonian ([T]) with the aid of 
time-dependent coupled mode theory ^19j |. In particular, 
we decompose the rii region, which behaves as an open 
cavity, into a set of modes whose evolution (assuming 
strong confinement in the cavity) is found to be: 



dak 
dt 



-c{t), fee (4) 



being the amplitude of the k-th cavity mode, with 
energy la^p, resonant frequency ujk and total lifetime 
= 7^ + 7- , the latter composed by the intrinsic decay 
rate l/rko and the additional escaping rate l/rg due to 
the coupling with the external photons source c{t). In 
order to model a reservoir of photons whose frequencies 
satisfy kL 1 we considered a broadband source c(t) = 
Q{t) J dwe'"*, which is switched on at t = 0. In this 
condition, the total electromagnetic energy H = |afc|^ 
stored into the cavity is found to be: 



n 



dco 



E 

k 



1 + e ^fc — 2 cos[(a;fc — LL>)t]e ^fc 



(5) 



Equation ([S|) allows to analytically investigate the har- 
vesting capacity of the system when the complexity C 
is varied. We begin when the photons dynamics at 
the short-wavelength scale shows the maximum unpre- 
dictability, i.e., for C = 1. In this condition, the long 
time dynamics [tp{t),9{t)] of photons will be distributed 
according to the natural measure [13, of the chaotic 
billiard, regardless of the photons initial position in the 
phase space. As a consequence, the rate of energy leakage 
will be independent of [■(/'(O), 0(0)], too. Energy leaks, in 
fact, relies on the transmission at the boundary S that, 
in turn, depends on the time-distribution of [^{t),6{t)] 
whose dynamics always converges to the same proba- 
bility function (i.e, the natural measure). The short- 
wavelength limit is an abstract scale that does not rely 
on a particular photons wavelength but equally applies 
to all wavelengths A (and frequencies w — 27rc/A) satis- 
fying ki ^1. As a result, when considering the more 
refined wave-envelope time-scale, we can conclude that 
under the maximum complexity degree the decay rates 
of the cavity modes should all be independent on the 
photon's frequency, thus implying 1/t/j = 1/r. In ad- 
dition to that, a second observation can be drawn from 
the Bohigas-Giannoni-Schmit conjecture Q, which states 
that the energy spectrum of a generic quantum (or wave) 
system with underling classical chaotic dynamics is dis- 
tributed in the same way as spectra of the standard ran- 
dom matrix ensembles within the same symmetry class. 
The chaotic billiard dynamics arising from the HJ equa- 
tions, in particular, belongs to the Gaussian Orthogonal 
Ensembles (GOE) class and should display repulsion in 
its wave spectrum, which will be characterized by non 
overlapping resonant frequencies ujk- Under this condi- 
tion, Eq. ([S|) can be further simplified as the integral 



yields significant contribution only for uj k, ujk'- 
T [i-e 

k 



nTe-Y.^\l-e-if ^UT^{l-e-i)\ (6) 



where the sum runs for all the M resonances Uk over- 
lapping with the spectrum of c(t). Equation ([5]), in its 
remarkable simplicity, predicts that under complex con- 
ditions all cavity resonances Uk cooperate in the accumu- 
lation of energy and sustain a coherent buildup process. 
At the steady state i ^> 00, Eq. © yields a form of the 
thermodynamic equipartition principle: 



dU r2 
— — = const — — 



(J) 



which predicts a linear energy increase in the "phase- 
space" given by the spectral domain w. Conversely, when 
the photons motion is completely predictable and the sys- 
tem complexity is zero, no natural measure exists for the 
whole phase space and different input conditions evolve 
toward different long-time probability distributions. In 
this regime, each mode resonance exhibits a different 
characteristic decay rate I/t^, and no cooperative dy- 
namics occurs. When the complexity falls outside these 
two limiting conditions, the phase space will be charac- 
terized by both predictable and unpredictable dynamics. 
In this situation, we expect Eq. ([B]) to be valid only in 
chaotic regions where a natural measure can be found, 
with the net result that only a small subset of mode res- 
onances can contribute to the energy buildup process. 
The impact of photons cooperation on the process of en- 
ergy harvesting can be evaluated from the power V that 
is transferred into the cavity. From Eqs. Q, this quan- 



tity is found to be P = 2 (1+^7^- 

The power V reaches the maximum when intrinsic and 
escape decay rates are equal, i.e., for rfco — t^, and 
rapidly goes to zero outside this limit. In the absence 
of any cooperative dynamics between photons, only a 
small number of resonances LOk could satisfy this condi- 
tion, while all the others will inevitably lead to a poor 
power rate transfer. As a consequence of that, the pro- 
cess energy storage is expected to show a saturation 
in the efficiency when many frequencies interact with 
the cavity, due to the limited amount of power that 
can be transferred under mismatched conditions. In the 
high complexity case, conversely, thanks to a frequency- 
independent mechanism of energy buildup, energy har- 
vesting is predicted to grow linearly with the source band- 
width thereby overcoming any bottleneck problem due to 
a mismatched load. 

Wave-carrier limit: ab-initio photons dynamics. — 
At the finest scale, we perform ab-initio simulations 
by numerically solving Maxwell's equations with paral- 
lel Finite-Difference Time-Domain (FDTD) simulations 
[2H . Ab-initio means "from first principles" and yields 
the most rigorous approach as it is based on the direct 
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FIG. 2. (Color Online) Ab-miUo results: (a)-(b) normalized 
energy H evolution for (a) single frequency excitation with 
A = Ao = 0.7 nm and (b) dense frequency comb spanning 
the range AA € [0.2, 1.2]^m; (c) energy maxima Hm versus 
source bandwidth A A (symmetrically increased from Ao) and 
complexity C (marker: FDTD simulation, solid line: numer- 
ical fit); (d) log- plot of the energy relaxation dynamics after 
switching off the source used in (b). 



solution of Maxwell's equations without any approxima- 
tion. Besides that, we also design a possible energy- 
harvesting device with application in the renewable en- 
ergy industry. Our device is constituted by a Silicon bil- 
liard in air {rio = 1), whose shape follows Eq. ([3]). At 
the input, we launch a dense frequency comb with wave- 
lengths A in the range AA € [0.2, 1.2]/Ltm, thus modehng 
a broadband source such as a supercontinuum or sun- 
light, and calculate the electromagnetic energy ^(C,^) 
stored into the cavity by a numerical computation of Eq. 
([1]) in the region delimited by rii (Fig. [2]). For a single 
frequency excitation with A = Ao = 0.7/im, no apprecia- 
ble distinction is observed as the complexity C is varied 
(Fig. ^jp^). However, as soon as the bandwidth AA of the 
input source is symmetrically enlarged from the central 
wavelength Aq, we observe a radically different evolution 
depending on the complexity of the classical motion of 
photons (Fig. (Jb-c). In the complex case C — 1, in 
particular, photons cooperation settles in and each cav- 
ity resonance ujk coherently contributes to harvest light 
energy in the system, thus leading to a linear increase 
of H. Conversely, when the classical motion of photons 
is characterized by a zero entropy, saturation occurs be- 
cause of the mismatch between intrinsic 1 /rfco and escape 
I/te decay rates of excited cavity modes (Fig. [Hd-c). To 
further investigate this aspect, we calculate the energy 
decay in the case of Fig. [2]d, i.e., when the input source 
spans the whole range AA G [0.2, 1.2]/iTO. This calcu- 
lation is performed by switching off the frequency comb 
after t = 16 ps, and then monitoring the time evolution 



of y, for At = 10 ps (Fig. [2ji). In perfect agreement 
with our predictions, a single exponential decay is ob- 
served when C — 1, while different decay times are 
manifested for C = 0. The effect of photons coopera- 
tion, as seen by comparing Figs [2^-b, is quite dramatic 
and leads to an energy increase of more than 65%. It 
is worthwhile remarking that such result is completely 
nontrivial, as the energy relaxation dynamics for C — 1 
is much faster than in the absence of complexity (see Fig. 
[2ji), and it can be predicted only by employing a multi- 
scale approach that keeps into account all scales of pho- 
tons dynamics. Besides that, it is not difficult to assess 
the impact of this result: nearly 50% of the price of any 
solar module lies in the cost of the Silicon; with such a 
technology we foresee new harvesting architectures were 
just a Si reshaping (with a constant volume) yield a dra- 
matic increase in the energy stored in the system, with 
a consequent breakdown of the price to produce a single 
energy Watt due to a larger energy production. 

Conclusions. — We have studied the problem of en- 
ergy harvesting in complex systems, by employing both 
analytic theory and ab-initio simulations. Time depen- 
dent coupled-mode theory and hamiltonian billiard dy- 
namics predict the existence of a mechanism of photons 
cooperation, which is entirely sustained by complexity. 
This phenomenon, stemming from the existence of a nat- 
ural measure in the billiard phase-space and verified by 
first-principle simulations, is observed as a dramatic in- 
crease in the electromagnetic energy stored into the sys- 
tem. In complete analogy to classical statistical mechan- 
ics where microscopic chaos predicts a single direction for 
the entropy, unpredictability in photons dynamics yields 
a single direction for the energy, which grows with the 
complexity degree of photons classical motion. This re- 
sult is expected to stimulate further theory and new ex- 
periments devoted to the realization of novel, complexity- 
driven, light harvesting schemes. 

We acknowledge discussions with G. Ruocco and P. de 
Bernardis. 



* landrea.fr atalocc hi@kaust.edu. sal | www.primalight.org] 
[1] D. T. Haar, Rev. Mod. Phys. 27, 289 (1955)| 
[2] J. Loschmidt, Sitzungsber. der kais. Akad. d. W. Math. 
Naturw. II 73, 128 (1876). 
P. Gaspard et al, Nature 394, 865 (1998). 
J. Halliwell, J. P. Mercader, and W. H. Zurek, eds.. Phys- 
ical Origin of Time Asymmetry (Cambridge University 
Press, Cambridge, 1994). 

S. W. Hawking, [Phys Rev. D 32, 2489 ( 1985) 

F. Haake, Quantum Signatures of Chaos (Springer, New 
York, 2001). 

G. Boffetta et al, Physics Reports 356, 367 (2002), 
G. Gallavotti, Statistical Mechanics: A Short Treatise 



[3] 
[4] 



[5] 
[6] 



[7] 
[8] 

(Springer, New York, 1999). 
[9] E. Gerstner, Nature Physics 4, 166 (2008). 



5 



[10] B. E. White, Nature Nanotechnology 3, 71 (2008). 

[11] P. J. Turnbaugh et ai, Nature 444, 1027 (2006). 

[12] J. L. Herek et ai, Nature 417, 533 (2002). 

[13] F. Caycedo-Soler et al, 

| Phys. Rev. Lett. 10 4, 158302 (2010) | 
[14] F. Cottone, H. Vocca, and L. Gammaitoni, 

Phys. Rev. Lett. 102, 080601 (2009) 
[15] P. Erdi, Complexity Explained (Springer, New York, 

2007). 



[16] E. Ott, Chaos in Dynamical Systems (Cambridge Uni- 
versity Press, Cambridge, 1997). 

[17] N. Chernov and R. Markarian, Chaotic Billiards (Amer- 
ican Mathematical Society, Providence, 2006). 

[18] S.-Y. Lee et ai, Phys. Rev. Lett. 93, 164102 (2004) | 

[19] H. A. Haus, Waves and Fields in Optoelectronics 
(Prentice-Hall, Engiewood Cliffs, N. J., 1984). 

[20] C. Beck and F. Schlogl, Thermodynamics of chaotic sys- 
tems (Cambridge University Press, Cambridge, 1993). 

[21] A. Taflove and S. C. Hagness, Computational Electro- 
dynamics: the finite-difference time-domain method, 3rd 
ed. (Artech House, 2000). 



